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We propose a microscopic model of critical current noise in Josephson-junctions based on individ- 
ual trapping-centers in the tunnel barrier hybridized with electrons in the superconducting leads. 
We calculate the noise exactly in the limit of no on-site Coulomb repulsion. Our result reveals a 
noise spectrum that is dramatically different from the usual Lorentzian assumed in simple models. 
We show that the noise is dominated by sharp subgap resonances associated to the formation of pairs 
of Andreev bound states, thus providing a possible explanation for the spurious two-level systems 
(microresonators) observed in Josephson junction qubits [R.W. Simmonds et al, Phys. Rev. Lett. 
93, 077003 (2004)]. Another implication of our model is that each trapping-center will contribute 
a sharp dielectric resonance only in the superconducting phase, providing an effective way to val- 
idate our results experimentally. We derive an effective Hamiltonian for a qubit interacting with 
Andreev bound states, establishing a direct connection between phenomenological models and the 
microscopic parameters of a Fermionic bath. 

PACS numbers: 74.50.-|-r 74.40.-|-k 



I. INTRODUCTION 

The performance of Josephson-junction devices func- 
tioning as units of quantum memory or as qubits depends 
to a large extent on the amount of charge and criti- 
cal current noise affecting each Josephson-junctiouFii^i^ 
One mechanism for critical current noise is to assume 
that trapping-centers (TCs) located in the tunnel bar- 
rier will partially block conduction whenever they cap- 
ture electrons from one of the superconducting electrodes 
[Fig. [Ua)]i^ The noise resulting from each TC is tra- 
ditionally modeled as two-level telegraph noise, with a 
Lorentzian noise spectrum, and a combination of several 
TCs leads to 1/f noisci^ 

Nevertheless, sensitive spectroscopy experiments on 
current-biased Josephson-junctions (phase qubits) re- 
vealed the presence of a few microwave resonators on top 
of the expected 1/f noisci^ These microresonators behave 
as spurious two-level systems buried within the tunnel 
barrier, whose coupling to the qubit produces reduced 
measurement fidelity and decoherencei^ Similar microres- 
onators were observed in flux qubits The microscopic 
origin of the microresonator remains unknown. However, 
there is strong evidence that improving the junction ox- 
ide quality reduces their concentration^ Recently, phe- 
nomenological models based on resonant coupling with 
the Josephson energy;^ and dielectric resonancei^ were 
proposed. Two measurement schemes to distinguish 
these different models were suggested J^ii^ An interest- 
ing connection between the low and high frequency scales 
of the noise spectrum due to a large number of microres- 
onators was demonstrated!^ To our knowledge, there are 
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FIG. 1: (a) A current biased Josephson junction adversely 
affected by the tunneling of electrons between one of the su- 
perconducting leads (sc) and a trapping-center defect in the 
insulating barrier (ins.). Josephson junction critical current 
noise is directly related to fluctuations in trapping-center oc- 
cupation due to modulation of the tunneling rate between su- 
perconducting leads, (b) A Cooper-pair box quantum bit af- 
fected by charge noise produced by a single trapping-center in 
the barrier, (c) Proposed set-up for measuring trap noise close 
to the superconducting transition. A single electron transistor 
(SET) is weakly coupled to an artificial trap, e.g. a normal 
state quantum dot or a nanotube. 



two proposals in the literature for the microscopic origin 
of these microresonators. The first is based on macro- 
scopic resonant tunneling in large Josephson-junctionSf^ 
This model explains the splitting of the Josephson en- 
ergy but predicts no dielectric resonance for the mi- 
croresonator. The second microscopic model is based 
on the structural two-level system in glassesjiSii^ This 
gives rise to the same dielectric resonance above and be- 
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low the superconducting critical temperature. Recently, 
a quantum computer architecture using microresonators 
as qubits was proposed.— Therefore, understanding the 
microscopic origin of the microresonator is of central im- 
portance for improving superconducting qubits. 

In a previous letter, we studied the charge noise 
spectrum due to a single TC hybridized with a non- 
superconducting Fermi sea.— At high temperatures we 
showed that the presence of a single TC with energy 
level close to the Fermi level leads to the expected 
Lorentzian spectrum characteristic of semiclassical ran- 
dom telegraph noise. At lower temperatures and frequen- 
cies below the TC linewidth, the noise has a quantum 
Johnson-Nyquist form reflecting the electron-hole exci- 
tations in the gate electrode Fermi sea.— 

Here we consider the noise spectrum due to a single 
TC hybridized with electrons in a superconductor. We 
show that the noise spectrum of each TC is characterized 
by a sharp resonance associated to the Andreev bound 
states formed from the TC hybridization with the super- 
conductor. We further demonstrate that our theory de- 
scribes a direct connection between this TC physics and 
the spurious microresonators observed at subgap frequen- 
cies in Josephson-junction devices.— i^i^ We derive an ef- 
fective Hamiltonian connecting the discrete levels to the 
microscopic parameters of the TC plus superconductor 
bath. Since the TC has an electric dipole moment, our 
model predicts that a sharp dielectric resonance will ap- 
pear when the sample becomes superconducting. 

The model proposed here is based on tunneling events 
between individual TCs and the superconductor. This is 
different from other model a^^'^^i^° that considered charge 
tunneling between two TCs mediated by Andreev states, 
resulting in a smooth noise spectrum that does not give 
rise to microresonators. 

Trapping-center fluctuation nearby to single electron 
tunneling devices^! is also an important source of charge 
noise and decoherence of charge qubits such as the 
Cooper-pair bo x^^i^^ [Fig. [T]Jb)] and double quantum 
dot.^ 

Fig- IHc) suggests a test device to probe TC noise 
around the superconducting transition temperature Tc 
that allows verification of our predictions in a controlled 
manner. A tunable artificial trap, which can be realized 
by a quantum dot or a nanotube in the normal state,— 
is coupled to a large metallic reservoir at temperatures 
close to the superconducting transition. A single elec- 
tron transistor (SET) is proposed to measure TC charge 
occupation in real timej^ This will map the emergence 
of the subgap resonance as T is lowered below the super- 
conducting transition temperature Te- 



ll. QUBIT DECOHERENCE AND QUANTUM 
NOISE 



Hence one can design circuits that behave as artificial 
two-level systems, realizing promising qubits for scalable 
quantum computer architectures.—. 

Consider a model for an artificial two-level system. 
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"^Qubit — 77^^(-^c) • 



(1) 



where & — {cTx,ay,az) is the vector of Pauli matrices 
denoting the qubit, and is a vector with dimensions 
of frequency. The latter is a function of Ic, the critical 
current of one of the Josephson-junctions in the circuit. 
We assume the critical current depends on the quantum 
state of TCs in the barrier, hence we write it as an opera- 
tor (notation Ic to distinguish quantum operators from c- 
numbers such as Ic). For small fluctuations we may write 
fliic) w flo + n'o{Sic)+0{5ic)'^, where flo fti{Ic)) and 
((5/c) = Ic — {Ic)- Choosing a coordinate system with z- 
axis along Hq, and a;-axis along fli = fig — (i7g • z)z we 
get 



^Qubit = -^fi^o^z + hr]:z{SIc)az + hri^{5Ic)crx- (2) 

Fluctuations in 5Ic affects the qubit through the param- 
eters r]z = ^J^o ' ^ ^'^^ ~ 5^1- "Fhe former leads to 
phase relaxation or decoherence, while the latter causes 
energy relaxation. 

In the weak coupling regime, all relaxation effects are 
fully characterized by the critical current noise spectrum. 



Si{io) = 



dt 



ic{t)-{ic) ic{o)-{ic) 



(3) 

For example, if the qubit is prepared in the excited state 
(It)), its rate of approach towards thermal equilibrium 
will be given by 
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(4) 



Similarly, if the qubit is prepared in a superposition state 
(IT)+ its coherence envelope = {g^ + i<^y) 

will be affected by low frequency noise according to ^^'^^ 



|(a+(t))Hcxp 



duj Si{uj)T{t, u)) 



(5) 



Here the filter function J-'{t, lu) depends on the particu- 
lar method chosen for probing qubit coherence. For free 
induction decay we have 



1 sm^ 
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(6) 



while for the Hahn echo 



The behavior of superconducting circuits containing 
Josephson-junctions is markedly quantum-mechanical. 
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with qubit coherence probed a.t t = 2te after the appHca- 
tion of a 7r-pulse at time tg.- Note that ^Haiin(2ie, 0) = 0. 
The Hahn echo filters out terms proportional to 5/(0) 
hence leading to much longer coherence times for qubits 
subject to low frequency noise (See Ref. [2^ for further 
discussion and derivations). 

For the purposes of this work, it is instructive to use 
Eq. ([5]) to study the effect of a sharp frequency peak (a 
resonance) in the noise spectrum. Assume Si{uj) has a 
sharp peak centered at Qucs with linewidth 1/tc;, 



Td 



1 



(8) 



Using Eqs. ([5]) and ([6]) and assuming ilRcs ^ l/Td we get 



exp 



Res 



cos fiRos^ 



(9) 

Therefore a resonance in the noise spectrum leads to loss 
of visibility of coherence oscillations. The loss of visibility 
is initially oscillatory, but decays exponentially to a fixed 
contrast for t ^ Td, similar to [30|. Although Eq. Q 
was calculated for free induction decay, it is also a good 
approximation for Hahn echoes in the limit JIros 3> 1/ie- 
The above discussion makes clear the fact that the key 
quantity to be studied in the context of qubit relaxation 
and decoherence is the time ordered noise spectrum de- 
fined by Eq. ([3]). If noise is the object of interest, the 
qubit acts as a spectrometer for quantum noisci^ Later, 
in section IVIII we are going to show that the same basic 
Hamiltonian also leads to the formation of avoided cross- 
ing with Andreev levels acting as junction resonators. 



III. MICROSCOPIC MODEL FOR CRITICAL 
CURRENT NOISE 

A. Trapping-center model Hamiltonian 

The Hamiltonian for a trapping-center coupled to a 
lead with Bardeen-Cooper-Schrieffer (BCS) interactions 
is given by^^i^ 



H ^Hq + Hbcs + V. 
The unperturbed trap Hamiltonian reads 



(10) 



(11) 



where rio- = d'ldcr is the electron number operator for a 
TC with spin tr i: ^^nd dj. is a Fermion creation oper- 
ator. The trap energy level Cd is measured with respect 
to the Fermi level (we assume ep = 0)- The unperturbed 
mean-field Hamiltonian for a superconducting Fermi lead 
is given by 

Wbcs = J2 ''k^ka - J2 ^cltctfci + H-C-, (12) 

k,a k 



where cj.^. creates a conduction electron in the gate elec- 
trode with energy e^, and Uka = c\^Cka- A is the super- 
conducting order parameter. The conduction electrons 
are hybridized with the TC via the hopping Hamiltonian 



(13) 
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where Vk is the tunneling matrix element for the electron 
between the TC and the superconducting lead. 

Here we assume TCs for which the on-site Coulomb re- 
pulsion, of the form t/nfn|, can be neglected. We remark 
that the chemical structure of TCs in the Josephson bar- 
rier is not known. There are many possible kinds of TCs 
associated with the amorphous oxide in a typical Joseph- 
son junction: 0-H complexes, various kinds of vacancies, 
dangling bonds, etc. Our model will be applicable to TCs 
with [/ ^ A. The U = Q idealization is an important 
starting point, particularly because it allows an exact so- 
lution of the noise problem. As we show below, our model 
seems to explain some of the important features observed 
in spectroscopy of Josephson qubits. In section rVIIII we 
discuss the expected modifications when U > {). 



Trapping-center fluctuation as a mechanism for 
critical current noise 



We now describe a model for the effect of TC fluc- 
tuation on the critical current of a Josephson-junction 
[Fig. [Ija)]. Our aim is to establish a direct relationship 
between critical current noise and the TC noise spectrum: 



S'„(cj) = 



|^e-*([n(t)-2n]K0)-2n]). (14) 



Here n = d^d^ is the total number operator for 
electrons occupying the TC level. The notation {A) = 
Tr {pgA} denotes grand-canonical averages using the 
density operator pc = e^'^'^^^^V^G, with (3 = l/fcsT 
and Zq the grand canonical partition function. In the 
absence of a magnetic field, = = n, therefore 
we write {n) = 2n to simplify the notation. 

The presence of a TC will produce weak modulations 
on the junction potential barrier F^i^ii^ Our model is to 
assume that the channel average matrix element for elec- 
trons tunneling from one lead to the other depends on n 
according to 



Tlr ~ T, 
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(15) 



The critical current Ic [or equivalently, the Josephson 
energy Ej = {h/2e)Ic\ is proportional to the modulus 
squared of Eq. (jlSP i^ so that in the adiabatic limit, 
for frequencies smaller than the inverse tunneling time^^ 
this directly translates into a fluctuation of the critical 
current 
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The critical current noise is therefore given by 

dt 



7-2 |7-(1)|2 
\2 



2ti 



e"^*([n(t)-2n][n(0)-2n]) 
(17) 



Hence within the linear approximation [Eq. (|15p ] the 
resulting critical current noise is directly proportional to 
the TC charge noise, S'„([j). The proportionality con- 
stant can be extracted directly from experiments probing 
critical current noisCi ^'^'^^ Below we focus on theoretical 
calculations of the TC noise spectrum Snioj) under dif- 
ferent parameter regimes. 

Our model assumes the TC is coupled to only one of 
the superconducting leads. Within the one-lead approx- 
imation critical current modulations are assumed to oc- 
cur only through variations of inter-lead tunneling due 
to population/depopulation of the trap [Eq. US])]. We 
therefore neglect the possibility for the trap electron to 
enter through one lead and exit through the other. These 
processes will lead to interesting phase dependent effects 
in the Josephson current.— We are not aware of stud- 
ies of critical current noise in this regime. Nevertheless, 
for zero phase difference between the leads, we may map 
the problem into a TC coupled to a single lead.— There- 
fore our results should remain valid in this case provided 
the phase is set to zero. In an experimental sample con- 
taining a few TCs we should expect that some of these 
are coupled to a single lead, others are coupled to both 
leads. The former case will lead to phase-independent 
noise, while the latter is expected to generate a phase 
dependent noise spectrum. In this context the theory 
developed here should be compared to measurements of 
the phase independent contributions to critical current 
noised Note that the TC only couples to both leads if 
it is in the middle of the junction with a difference in 
separations to either lead being smaller than a tunnel 
length. Given that junctions are typically much thicker 
than a tunneling length to the extent that the latter are 
knownf2£ the present theory covers most of the possible 
TC locations. 

In this work we calculate the noise spectrum under the 
assumption that the TC remains in thermal equilibrium 
with the superconducting reservoir. Therefore our results 
are valid at the regime where non-equilibrium effects are 
weak or can be neglected. This is the case for a current- 
biased Josephson-junction in the zero voltage state, or 
whenever the voltage is low enough so that the electrons 
in the lead may still be characterized by a Fermi distribu- 
tion. The thermal equilibrium assumption implies that 
the noise spectrum satisfies the detailed balance condi- 

tion, S{—uj) = e 3(^10). The finite frequency noise 
spectrum measured by a particular detector depends on 
details such as the detector temperature To (not neces- 
sarily equal to the TC plus Fermi sea temperature T). 
For example, current noise measured by an LC circuit 
relates to our calculated time ordered noise [Eq. (fTT)) ] in 



the following way^ 



(18) 

where K denotes the effective coupling constant between 
the current carrying wire and the LC circuit. The exper- 
iment proposed in Fig. [TJc) should be interpreted using 
Eq. (HID. 



IV. RELATIONSHIP BETWEEN NOISE AND 
TRAPPING-CENTER SPECTRAL FUNCTIONS 

In this section we show that the TC noise spectrum 
Eq. can be expressed as an integral over all possible 
quasiparticle-quasihole excitations in the TC plus super- 
conductor problem. In order to derive this result, we 
define the Matsubara and real time correlation functions 
as follows^ 

5(t-t') = -Tr{pGT,[<5?i(T)a(T')]}, (19a) 
S'(^)(t-t') = ~ie{t~t')TT{pG[5nit),5n{t')]},{19h) 

where we used the notation 6n = n — 2n. Here, we 
use the Matsubara representation of operators n^r) = 
gTW/^^^Q-jg-rW/;!^ that are obtained from the Heisenberg 
representation by substituting it r. 

Applying Wick's theorem to Eq. (|19ap leads to 

(20) 

where we have introduced the normal Q and anomalous 
T TC Matsubara Green's functions, 

Q„„.[t) = -Tr{/5Gf,[d,(T)4,(0)]}, (21a) 
T„„.{t) = -Tr{/5GTr[da(T)d.'(0)]}. (21b) 

We now take the Fourier transform of Eq. (|20p . 

S{iujn) = ^11^ dre^'^'^'^ S(t), and insert the Lehmann rep- 
resentation for the TC Green's functions, 



T„„,{iuJn) = h 



dio' 
duj 



lOJn — ^ 



(22a) 
(22b) 



The TC spectral functions Aaa' {uj) and Baa' {^) play a 
fundamental role in our theory. For a BCS model such 
as Eq. pU]) . we have A]] = An = A with A real, and 
A]i ~ An ~ 0. Also, the spectral function related to 
Gorkov's T function is non-zero only for = = B, 
with B real. After inserting these Lehmann representa- 
tions into Eq. (PO)) . the result is readily evaluated using 
the residue theorem, and taking advantage of the fact 
that S{iujn) is non-zero only at even (Bose) Matsubara 
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FIG. 2: Particle-like Andreev bound state as a function of 
hybridization 7 for different trap energies td, all in units of 
the superconducting gap A. The hole-like Andreev bound 
state has the same energy but opposite sign. 



Eq. (j25bp shows that the spectral function ,S(e') describes 
the extent to which TC pairing is induced through its 
hybridization with the Fermi gas, i.e., S(e') can be in- 
terpreted as a single-state proximity effect. Integrating 
Eq. 123]) over all frequencies and using Eqs. (|25ap . (j25bp 
we obtain the following sum rule, 



Sn{uj)duJ 



E 



(26) 

As a cross check, Eq. can be derived directly without 
using Eq. (|23|) by simple applications of Wick's theorem 
and use of the Fermi anticommutation relation, such that 
= rio-. Interestingly, the sum rule [Eq. ([26|l ] shows 
that the onset of superconductivity tends to increase the 
amount of noise produced by a TC. 



frequencies [cun = wk / {h(3) with n even]. Finally, analytic 
continuation {iujn — > w -|- irf) allows us to extract the TC 
noise spectrum from the imaginary part of S'^^^(cj). This 
leads to a convenient expression for the TC noise spec- 
trum, 

/oo 
de' [A„,,{e')A,,,{e' -uo) 

~B;,Ae')B.'A^' - ^)] [1 - .fie')]fie' - c.)(.23) 



Here the Fermi functions are given by 

1 



(24) 



The expression Eq. (^5]) is an exact result. Its deriva- 
tion relied on the use of Wick's theorem, which is valid 
only for a quadratic Hamiltonian Eq. PU)) (C7 = 0)4i 
It expresses the fact that the TC noise spectrum is the 
sum of all quasiparticle-quasihole excitations involving 
the dressed TC plus superconductor at thermal equilib- 
rium. Eq. (j23p is the generalization of an equation de- 
rived by us previously, using a canonical transformation 
in the TC plus normal metal Fermi sea problem (See 
Eq. (8) in Ref. [H). 

It is instructive to derive a sum rule for the noise spec- 
trum starting from Eq. (|23|) . First, note that the TC oc- 
cupation number and TC pairing correlator are related 
to the spectral functions in the following way, 

= {dUa)= j de'A{e')f{e'), (25a) 

Fa = -(dT^i) = / de'B{e')f{e'). (25b) 



V. TRAPPING-CENTER SPECTRAL 
DENSITIES AND ANDREEV BOUND STATES 



The TC spectral densities are known exactly for the 
case of zero on-site Coulomb repulsion.— 1^1^ These can 
be written as 



A = Ae{\e\ - A) + [a+5{e 
B = B6{\e\- ^)+b+[5{e-Eb) 



Ei.)+a^6{€ + Ei,)] ,(27a) 
Eb)-5{e + Eb)], (27b) 



where 6 is the step function. Each spectral density is 
composed of a continuous above-gap component which 
is non-zero only at energies outside the superconducting 
gap (|e| > A). For energies within the gap, there are 
two sharp Andreev bound states, with positive particle- 
like energy E}, and negative hole-like energy —Eb- These 
bound states are reminiscent of the TC localized level e^, 
whose energy is renormalized to ±£^f, due to hybridization 
with Cooper pairs. 

In order to express these functions analytically, we de- 
fine the trap hybridization parameter 7 as 



7 = ■ngo{V^)e 



(28) 



where Vk is averaged over tk 
density at the Fermi level. 



e_F and is the energy 
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The above-gap contributions are given by 
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2. Subgap resonance 
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FIG. 3: (Color online). Depiction of the energy density as 
a function of energy for the TC plus superconductor model, 
and the most important quasiparticle-quasihole excitations 
(denoted by arrows) determining the TC noise spectrum. 



Note that 0+ 7^ a_ in the asymmetric case ed 7^ 0, but 
6-1- = h- always. A useful relation is that 6+ = — ^a+a_. 
In Fig. [2] we plot the Andreev levels Ei, as a function of 
TC hybridization for different TC energies Cd- 



It is useful to establish a connection to the case of 
a point contact between superconductors4^ In this case 
the transmission of electrons across the point contact is 
dominated by the presence of two Andreev bound states 
at equal and opposite energies with respect to the Fermi 
level. For zero phase difference these Andreev levels are 
located close to ±A. In our case, the trapping-center is 
equivalent to a point contact provided je^l 3> A; Looking 
at Fig. [5] we see that is indeed slightly below A in this 
limit. 



where sgn (e) denotes the sign of e. We remark that 
these are finite temperature spectral densities; the tem- 
perature does not appear explicitly because the Matsub- 
ara Green's functions for [/ = depend on temperature 
only through the Matsubara frequencies4^ The Andreev 
bound state energy is given by the single pair of real roots 
±£^6 of 



£2 ( 1 



27 



VA2 - E"^ 



4-1^^ 0, 



(30) 



with the amplitudes a± and 6+ of Eqs. (|27a[) and (j27bp 
given by 



(A2 



El) 



a± = 



{cd ± Ebf + 7' 



2[(2A2-£;2) (e2 + ^2) 
-7Agfcv/A2-£;2 

[(2A2-£;2)(,2+^2)_£;4]- 



(31a) 
(31b) 



VI. EXPLICIT RESULTS FOR THE NOISE 
SPECTRUM 



We now show explicit results for the noise spectrum of 
a single TC hybridized with a superconductor. The an- 
alytic expressions for the spectral functions are inserted 
into Eq. where in the absence of a magnetic field 

c' = 2 and Aa-a = Baa = 0. For oj > the noise is 
given by 
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+6i(tj-2A) / de 



-oo jA+ij JA 

+2a+a_[l - fiEb)]f{-Eb)6iiJ - 2Eb) 



[Aie)A{e - lo) - B{e)B{e - u)] [1 - /(e)] /(e - a;) (32a) 

(32b) 



- ^Eb - A) [a^Ai^Eh a.^Ai-^Eb + u;) - 2h+B{£,Eb - uo)] [1 - f {£.Ei,)]J {£.Eb - uj) \ .(32c) 
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FIG. 4: (Color online). Trapping-center noise spectrum near 
the superconducting transition temperature Tc, for = 0. 
For T > Tc the noise has the Lorentzian form characteristic 
of random telegraph noise^^. As T is lowered below Tc a gap 
opens in the noise spectrum, and a sharp subgap resonance 
appears as a transition between two Andreev bound states 
(for convenience, we represent the subgap resonance with a 
phenomenological linewidth equal to O.OI7). This shows that 
TC noise is dramatically affected by superconductivity. 



The cij < expression can be obtained from detailed bal- 
ance S{-u}) = e~'^S{u;). 

The positive frequency spectrum is interpreted as the 
sum over all possible quasiparticle-quasihole pairs cre- 
ated vi^hen the TC plus Fermi sea absorbs a photon with 
energy huj emitted by the noise detector. Fig. [3] illus- 
trates the energy excitations associated to TC noise. The 
first contribution is a continuum-continuum transition 
Eq. p2ap where the hole (particle) is in the continuum 
below (above) the superconducting gap. This gives a 
smooth contribution to the noise spectrum when uj > 2A. 
The second line Eq. (j32bp is the subgap resonance. The 
resonance occurs when a hole is created at Andreev level 
~Eb and a particle is excited at level +Ei,. This contri- 
bution is a sharp transition between Andreev levels: The 
noise is a delta fimction peaked at a; = 2Eb. The third 
line Eq. (|32c[) refers to transitions involving one of the 
Andreev levels and the continuum. This gives smooth 
contributions for to > A ± Eb- 

Fig.|4]shows the noise spectrum for temperatures above 
and below the critical temperature for transition into the 
superconducting state. We assume ed ~ 0, with the su- 
perconducting energy gap dependent on temperature ac- 



cording to A 1.76kBTc^/T^^^T/T'c for T < T^, and 
A = for T > TcM. We assumed ksTc/j = 11.96, con- 
sistent with the value of = 1.196 K for aluminum with 
a trap hybridization parameter ^/kg = 0.1 K. For T > Tc 
the noise is a Lorentzian with linewidth 27 / h, consistent 
with the high temperature limit for random telegraph 
noise discussed in Ref. [l3l (Note that kgTc » 7 in 
Fig. 2]). As the temperature is lowered below Tc a sharp 
resonance appears at energy equal to two i^^.To display 
the subgap resonance in the figure we represented the 
delta function as a Lorentzian with linewidth equal to 
O.OI7. 

For T < Tc, kgT < Eb, and Cd < A the noise is well 
approximated by 

S'„(tj) « 2h{2a+a-5{uj - 2Eb) + e{uj - Eb - A) 
X [a+A{Eb -uj) + a-A{-Eb + uj) 
-2b+B{Eb-Lj)]}. (33) 

Fig.Oshows the low temperature noise spectrum (T <C 
Tc), with parameters normalized by the superconduct- 
ing energy gap. We also show the breakdown of the 
noise into its various contributions. For convenience, 
we represented the subgap resonance as a Lorentzian 
with linewidth O.OOlA. Our theory does not account for 
broadening mechanisms, but we expect that disorder and 
other inhomogeneities will be a source of broadening for 
Andreev bound states. 

For the parameters of Fig. [5] the subgap resonance 
accounts for 59% of the noise power. The remainder 
is due to Andreev-continuum transitions (33%), with 
continuum-continuum transitions contributing only 8%. 
Remarkably, the continuum-continuum contribution is 
quite small, in spite of being responsible for all the 
Lorentzian noise aXT > Tc (in the normal state). 

Fig. [6] shows the low temperature noise spectrum for 
a case where the Andreev bound states are very close to 
the gap edge, Eb = 0.981A (parameters = 0, 7 = IDA, 
and fcsT = O.IA). This case is quite different from Fig. (5] 
The continuum-continuum contribution is now 94% of 
the noise power, with Andreev-continuum contributions 
5.7%, and subgap resonance contributing only 0.3%. 

Fig. [7] depicts the noise in the asymmetric regime 
(ed ^ 0), with e<j = 5A, 7 = 0.5A, and fc^T = O.IA 
(the Andreev bound states are at Eb — 0.999A). Here 
the continuum-continuum contribution accounts for 98% 
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FIG. 5: (Color online). Trapping-center noise for tempera- 
tures well below the superconducting transition. The param- 
eters are 7 = 0.5A, = 0, fcsT = O.IA, leading to Andreev 
energy Et = 0.35A. About 60% of the noise power is due to 
one sharp subgap resonance represented here by a Lorentzian 
with linewidth O.OOIA. The remaining 40% is dominated by 
processes involving the creation of a hole in the continuum 
and the excitation of an Andreev level at +Et. This occurs 
only for o) > A -I- Et- 
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FIG. 7: (Color online). TC noise spectrum for the asymmet- 
ric case ed = 5A, with Andreev levels at ±Eb = ±0.999A. 
The noise spectra has a smooth peak close to hu — 6 A. This 
occurs because the spectral function peaks at ed- Similar to 
Fig. the noise is dominated by continuum-continuum con- 
tributions (98%) with the subgap resonance contributing less 
than 0.1%. 



VII. ANDREEV STATES AS JUNCTION 
RESONATORS 
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FIG. 6: (Color online). TC noise spectrum in the supercon- 
ducting regime, in a case where the Andreev bound states are 
very close to the gap edge {Et = 0.981 A). Here the noise 
is dominated by continuum-continuum contributions (94% of 
the noise power) , with the subgap resonance contributing only 
0.3%. 



of the noise, with Andreev-contimmni transitions con- 
tributing « 2% and subgap resonance contributing less 
than 0.1%. Interestingly, the noise has a broad peak at 
hio = 6 A, that occurs because the spectral functions have 
a smooth peak at — 5A. 

Figs.[M7]show that the noise changes its character com- 
pletely due to the opening of a gap and the formation of 
Andreev bound states in a superconductor. 



We now relate our theory to the experimental obser- 
vation of "spurious two-level systems" (microresonators) 
in phase-based^ and flux-based^ superconducting qubits. 

The model Hamiltonian for the interaction of a qubit 
with a TC plus Fermi sea is simplified by projecting 
onto the Hilbert space of Andreev bound states. This 
is achieved by expressing the TC operators as 



V2 



V2 



4 



V2 



V2 



ifT,cont, (34a) 

" rfl.cont, (34b) 



where alj. is a creation operator for an Andreev level with 
energy zLEb, and the da-^cont denote the additional opera- 
tors acting on the continuum. The canonical transforma- 
tion defined by Eqs. (|34ap and (j34b[) diagonalizes our TC 
model when u = Y^a+ and v = ^/aZ. This can be verified 
by calculating the Green's function using the canonical 
transformation and comparing to Eqs. ((27a|) and ([m]) . 
Substituting Eqs. ((Ma|) and (|34bp into Eq. ^ wc get 
an effective Hamiltonian for the qubit interacting with a 
pair of Andreev bound states, 

7Yq_a ^ ^fi^o^z + Eba\_a+ - Ei,a^_a_ 

(35) 



+ {a+ - a_) (a\a 



Here = h{5Ic)'riz and A^. — h{SIc)'>]x arc characteris- 
tic coupling energies between the qubit and the Andreev 
levels. For phase qubits these should be a fraction of the 
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change in Josephson energy h{5Ic)/{'2,e). Recall from sec- 
tion [IT] that the Ty's depend on qubit design, while [5Ic) 
is the characteristic shift in critical current due to a TC. 
A similar expression will hold for other kinds of qubits, 
for example, in a Cooper-pair box ^ PqPtc / is the 
electrostatic energy due to the interaction of the Qubit 's 
electric dipole moment pg and the TC (dipole moment 
PTC due to the image charge produced at the reservoir) 

The qubit-Andreev interaction is weighted by addi- 
tional factors accounting for the branching of the impu- 
rity spectral weight into different channels — not all of 
the impurity's noise goes into the Andreev channel. The 
first interaction, 2y/a+aZ{ai^a- + aLa+){Xzd'z + Xx^'x) 
produces admixture between qubit and Andreev lev- 
els, and leads to important anticrossings in qubit spec- 
trometry. The second interaction, [a^ — a-){a\^a+ — 

a'_a-){Xz(Tz + Xx<^x) only exists in the asymmetric case 
{cd 7^ 0). It enables the design of quantum gates through 
electrical manipulation of Andreev states. 

The Hamiltonian Eq. p5|) describes a four-level sys- 
tem, where the qubit energy levels are hybridized with 
the pair of Andreev states; it serves as a starting point 
to study non-equilibrium effects for a qubit coupled to 
Andreev excitations. Fig. [H^a) shows the energy levels 
Ei obtained after diagonalizing Eq. ((35|) for A^; = 0.2Eb 
and Xz = 0. Note the level anticrossing when H^Iq = 2Eb. 
Fig. [DJb) shows the two lowest energy transitions mea- 
sured by qubit spectroscopy, Ei — Eq and E2 — Eq. We 
remark the similarity of our Fig. [8{b) to the experimen- 
tal data in Fig. 2(a) of Ref. @. For these frequencies the 
qubit is highly mixed with the Andreev excitation. 

Therefore each pair of Andreev levels acts as a mi- 
croresonator, with frequency in the range 2Ei) G (0, 2A). 
The anticrossing behavior occurs only when the qubit is 
in resonance with a transition between Andreev levels, 
i.e., when the qubit frequency coincides with a subgap 
resonance in the noise spectrum. 



VIII. DISCUSSION AND CONCLUSION 

In summary, we calculated the noise spectrum due to 
individual trapping-centers (TCs) hybridized with a su- 
perconducting lead. We showed that the opening of a 
gap and the formation of Andreev bound states change 
the character of the noise completely. At T < Tc, the 
noise is substantially different from the usual Lorentzian 
spectra assumed in simple models. 

In many cases the noise is dominated by a subgap res- 
onance related to transitions between Andreev bound 
states at energies ±Ei, reminiscent of the localized TC 
states. At T <C Tc, the subgap resonance may account 
for over half of the noise power (See Fig. [5]). The remain- 
ing noise power occurs only at huj > A + Ei,, giving a 
smooth gapped spectrum related to the excitation of an 
Andreev level into the continuum. 

We assumed a TC model with zero on-site Coulomb 
repulsion. As a result, the noise can be expressed ex- 




FIG. 8: (Color online), (a) Energy level structure of the ef- 
fective Hamiltonian for a superconducting qubit interacting 
with a pair of Andreev bound states (one hole-like at energy 
—Eb and one particle-like at +Eb). The Andreev-qubit cou- 
pling is assumed A — 0.2Eb. Anticrossing behavior occurs 
when the qubit energy splitting matches the subgap reso- 
nance, huj = 2Et. (b) Energy differences between the ground 
state and the first and second excited states. The result is 
remarkably similar to spectroscopy measurements on phased 
and flux qubits.— 



actly as an integral over TC spectral densities, which 
are known analytically. This constitutes a limiting case 
which provides a fully characterized reference point. We 
now discuss the expected role of TC Coulomb repulsion. 
Spectral densities for U > were calculated using the 

"113. 



numerical rcnornialization group method in Refs. 
For U > Q, the energy Ef, of each Andreev level is shifted, 
but the number of Andreev levels remains the same (one 
hole-like and one particle-like per TC).— 

At [/ = 0, the total noise power [Eq. (^5)) ] is apprecia- 
ble only if ed lies within the interval [—A, A] (or within 
MaxIfc^T, 7} of this interval. For T > this result is 
equivalent to the one found in Ref. [l3|)- An interest- 
ing open question is whether this result will change for 
U >Q. 

We derived an effective Hamiltonian for a supercon- 
ducting qubit interacting with a TC, showing that the 
qubit sees the TC as two Andreev levels. Anticrossing 
occurs when the qubit frequency is in resonance with the 
energy separation of the two Andreev levels. This gives 
a microscopic explanation for the experimental obser- 
vation of microresonators coupled to Josephson-junction 
devices. Simmonds and collaborators observed anticross- 
ing behavior at a number of frequencies in the spec- 
troscopy of Josephson-junction phase qubitsi^ii Plourde 
and collaborators observed a similar effect in the spec- 
troscopy of flux qubits 4 Kim et al^ observed avoided 
level crossings in the spectra of a Cooper-pair box. 

Our work establishes a direct connection between TCs 
in the Josephson-junction insulator and the presence of 
these anticrossings. 

Another interesting implication of our model is that 
each TC will become a sharp dielectric resonance only 
when the lead becomes a superconductor. TCs are 
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charged defects, possessing an electric dipole moment due 
to their image charge in the superconducting lead. The 
fluctuation-dissipation theorem implies that the power 
absorbed by a TC irradiated by an AC electric field at 
frequency w is given by P^^ oc LuSn{i-o). Therefore, the 
subgap resonance in S',i((jj) can be detected as a sharp res- 
onance in dielectric absorption P^^ (in the normal state, 
Puj will be a broad resonance, see Fig. d]). 

This effect provides a powerful method to validate our 
theory experimentally. There are two other microscopic 
models for the microresonator: Macroscopic resonant 
tunnelingi^ results in no dielectric resonance; structural 
two-level systemiSiii^ gives rise to the same dielectric res- 
onance above and below Tc- Hence microwave absorp- 
tion experiments above and below Tc will clearly reveal 
whether the microresonator is a pair of Andreev levels or 
not. 

In conclusion, we have developed a microscopic the- 
ory for critical current and charge noise in superconduct- 
ing devices based on a charge tunneling model with indi- 
vidual trapping-centers. We showed that the supercon- 



ducting gap and the formation of Andreev levels plays a 
prominent role in determining the noise spectrum, pro- 
viding a microscopic explanation for the microresonators 
observed in experiments. Our calculated noise spectrum 
is drastically different from the usual phenomenological 
Lorentzian and 1/f noise spectra derived in previous 
work. 
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